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Abstract 

The gravitational interaction is scale-free in both Newtonian gravity and general theory of 
relativity. The concept of self-similarity arises from this nature. Self-similar solutions reproduce 
themselves as the scale changes. This property results in great simplification of the governing 
partial differential equations. In addition, some self-similar solutions can describe the asymptotic 
behaviors of more general solutions. Newtonian gravity contains only one dimensional constant, the 
gravitational constant, while the general relativity contains another dimensional constant, the speed 
of light, besides the gravitational constant. Due to this crucial difference, incomplete similarity can 
be more interesting in general relativity than in Newtonian gravity. Kinematic self-similarity has 
been defined and studied as an example of incomplete similarity in general relativity, in an effort 
to pursue a wider application of self-similarity in general relativity. We review the mathematical 
and physical aspects of kinematic self-similar solutions in general relativity. 
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1 Introduction 



Scale-invariance is one of the most fundamental characteristics of gravitational interaction in both 
Newtonian gravity and general relativity. This implies that if we consider appropriate matter fields, 
the governing partial differential equations are invariant under scale transformation. Due to this 
feature of the governing equations, there are self-similar solutions, which are invariant under the scale 
transformation. Self-similarity assumption enables us to simplify the governing equations. Self-similar 
solutions have a wide range of applications in astrophysics. See Jl] for a recent review of self-similar 
solutions in general relativity. See for self-similarity in more general contexts. 

When a theory has no characteristic scale, we can expect scale-invariance of the theory. In New- 
tonian gravity, the gravitational constant G, with dimension M _1 L 3 T~ 2 , is the only dimensional 
physical constant in the field equations, where M, L and T denote the dimensions of mass, length and 
time, respectively. It is impossible to construct a physical scale only from G. In general relativity, 
there exists another physical constant c, which is the speed of light, with dimension LT^ 1 . In spite of 
these two dimensional constants, no characteristic length scale can be constructed from these physical 
constants. However, due to the existence of these two dimensional constants, general relativity is qual- 
itatively different from Newtonian gravity with respect to scale invariance. If we consider quantum 
gravity, the Planck constant h appears, with dimension MI 2 T _1 , so that there exists a characteristic 
scale /pi = G l l 2 h}l 2 / 'c 3 / 2 , which is called the Planck length. Therefore, in the quantum theory of 
gravity, it is plausible that the scale invariance of the theory is broken down. Hereafter in this review 
we focus on Newtonian gravity and general relativity. We follow the sign conventions of [HS] for the 
metric, Riemann and Einstein tensors. 

2 Self-similarity in Newtonian gravity 

Since Newtonian gravity postulates an absolute system of space and time, we can directly apply 
the general formulation of self-similarity to this system A solution is called self-similar, if a 
dimensionless quantity Z(t, x) made of the solution is of the form 



where x and t are independent space and time coordinates, respectively, and a(t) is a function of t. 
This implies that the spatial distribution of the characteristics of motion remains similar to itself at 
all times during the motion. If the function a(t) is derived from dimensional considerations alone, 
i.e., if it is uniquely determined so that x/a{t) is dimensionless, the self-similarity is called complete 
similarity or similarity of the first kind pQ. In more general situations, the characteristic length or time 
scale may be constructed by the dimensional constants in the system. Then, the function a{t) cannot 
be uniquely determined from dimensional considerations alone. In such cases, self-similarity is called 
incomplete similarity or similarity of the second kind pp. For example, when we have the constant 
sound speed c s and no characteristic scale, then a(t) is uniquely determined as a(t) = c s t. In this case, 
the similarity is called complete. However, when we have a characteristic length scale I besides the 
sound speed c s , then a(t) = l l ~ a {c s t) a is possible and the constant a may not be determined from 
the governing equations. In this case, the similarity is called incomplete. The constant a may be 
determined by boundary conditions. It should be noted that the dimensional constant could appear 
not only from governing equations but also from boundary conditions. 

Here, we give two important examples of completely self-similar solutions in Newtonian self- 
gravitating fluid mechanics. The basic field equations for spherically symmetric hydrodynamics of 
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a self-gravitating ideal gas in Eulerian description are given by 

+ ( , 3) 

^r + ^ = ' (2 - 4) 

— = V ft (2.5) 

where p,v,M and G denote the mass density, radial velocity, total mass inside the radial coordinate 
r, and gravitational constant, respectively. 

2.1 Isothermal gas 

First we consider an isothermal gas as a gravitational source. Since the isothermal gas is a relevant 
description of cold molecular clouds in galaxies, self-similar solutions have been intensively studied in 
Newtonian gravity in modeling the star formation process [451 1611 1621 1411 176j . It has been revealed that 
self-similar solutions play important roles in the gravitational collapse of an isothermal gas [271 1721 ISU] . 
The stability of these self-similar solutions have been studied jSHl OH HUH EH] • A new insight has been 
obtained in this system in the context of critical behavior in gravitational collapse |471 139j . 

For an isothermal gas that obeys p = c^p, where c s is the constant speed of sound with dimension 
LT _1 , it is impossible to construct a characteristic scale from c s and G. We introduce the dimensionless 
self-similar coordinate 

_ Cst_ 
r ' 

for self-similar solutions. Then we also introduce the dimensionless functions U, P and m: 



(2.6) 



v(r,t) = -c s U(r,t), (2.7) 

KM) = (2.8) 

M(r,t) = 4t^t) (29) 

We assume that the above-defined functions U, P and m depend only on z. From this assumption, 
equations (|2.2|) - ()2.5|) become 

(zu + \y — z* 

zP[2 - P(zU + 1)] 

" {zU+ l ) 2_ z 2 » ^ LL > 

m = P(U + l/z), (2.12) 

-z 2 m = P, (2.13) 

where the prime denotes the derivation with respect to z. The self-similar solutions for an isothermal 
gas are obtained from these ordinary differential equations. Self-similar solutions scale for the scale 
transformations i = at, r = ar as 

v (r,t) = v(r,t), (2.14) 

p(r,t) = ^M, (2.15) 

M(r,t) = aM(r,t), (2.16) 
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where a is a constant. The basic equations for self-similar solutions are singular at the center and at 
the point at which (zU + l) 2 — z 2 = is satisfied, which is called a sonic point. 



2.2 Poly tropic gas 

Next, we consider a polytropic gas as a gravitational source. A polytropic gas obeys the equation 
of state p = Kp' 7 , where 7 is the dimensionless constant called the adiabatic exponent and K is a 
constant with dimension M 1 " 1 L 3 ~ t ~ 1 T~ 2 . As in the isothermal gas system, it is impossible to construct 
a characteristic scale only from G and K if 7 7^ 2. For the exceptional case, 7 = 2, the system has 
a characteristic length scale I = ^/KjG but even in this case the self-similar variable z is uniquely 
constructed. Then, complete similarity is applicable to this system. Self-similar solutions in this 
system have been studied [77HH3]- The stability of these solutions have been studied [3*31 l3*i| . 
For the polytropic case, we introduce the dimensionless self-similar coordinate 



K(-t) 



2-7 



(4vrG)(T- 1 )/ 2 r' 

Then we also introduce the dimensionless functions U, P and m: 



(2.17) 



v (r,t) = -(4vrG) (1 - 7)/2 v / ^(-t) 1 - 7 C/(r,t), (2.18) 

P{r,t > ~ (4vrG) 1 /(2-7) r 2/(2- 7 )' { - 2A9 > 
, K 3 / 2 (-t) 4 -^m(r,t) 

M ^ = (4 7 r)3(,-l)/2 g (37-l)/2 - ( 2 ' 2 °) 

We assume that the above-defined functions U, P and m depend only on z. In the polytropic case, 
the sonic point is defined by (2 — 7 — zU) 2 — ^z 2 ^ 2 ~^ = 0. Self-similar solutions scale for the scale 
transformations i = at, f = a 2 ~^r, as 

v(f,t) = a l ~^v{r,t), (2.21) 

M{r,t) = a 4 ~ 37 M(r,t), (2.23) 

where a is a constant. In this case, the scaling rates for r and t, which keep z constant, are different 
from each other. 

It should be again emphasized that in both the isothermal and polytropic cases, the self-similarity is 
complete since the self-similar variable z can be obtained from dimensional considerations alone. This 
is because there are only two dimensional constants in the system, while there are three independent 
dimensions M, L and T. 



3 Self-similarity in general relativity 
3.1 Homothety 

In general relativity, the concept of self-similarity is not so straightforward because general relativity 
has general covariance against coordinate transformation. This implies that the definition should be 
made covariantly in general relativity. In the following, we use units where the speed of light c is 
unity. In this choice of units, T = L is obtained and the velocity is dimensionless. 

In general relativity, the term self-similarity can be used in two ways. One is for the properties of 
spacetimes, the other is for the properties of matter fields. These are not equivalent in general. The 
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self-similarity in general relativity was denned for the first time by Cahill and Taub jllj . Self-similarity 
is denned by the existence of a homothetic vector £ in the spacetime, which satisfies 

= 2a 9fiu, (3.1) 

where g^ u is the metric tensor, denotes Lie differentiation along £ and a is a constant This is 
a special type of conformal Killing vectors. This self-similarity is called homothety. If a, ^ 0, then it 
can be set to be unity by a constant rescaling of If a = 0, i.e. C^g^ u = 0, then £ is a Killing vector. 

Homothety is a purely geometric property of spacetime so that the physical quantity does not 
necessarily exhibit self-similarity such as L^Z = dZ, where d is a constant and Z is, for example, the 
pressure, the energy density and so on. From equation (j3.1j) it follows that 

C^ uap = 0, (3.2) 

and hence 

£^iV = 0, (3.3) 

C^G^ = 0. (3.4) 

A vector field £ that satisfies equations (|3.'2[) . (|3.3f) and l)3.4[) is called a curvature collineation, a Ricci 
collineation and a matter collineation, respectively. It is noted that equations (j3.2j) . (|3.3j) and (|3.4|) 
do not necessarily mean that £ is a homothetic vector. We consider the Einstein equations 

G^ = Z-kGT^, (3.5) 

where T^ v is the energy-momentum tensor. If the spacetime is homothetic, the energy-momentum 
tensor of the matter fields must satisfy 

CgT^ = 0, (3.6) 

through equations ()3.5|) and ()3.4j) . For a perfect fluid case, the energy-momentum tensor takes the 
form of 

T^u = (p + ^)u^u v + pg^, (3.7) 

where p and fi are the pressure and the energy density, respectively. Then, equations 1)3. lj) and (|3.6|) 
result in 

= -au'\ (3.8) 

C^n = — 2a/x, (3.9) 

= -2ap. (3.10) 

As shown above, for a perfect fluid, the self-similarity of the spacetime and that of the physical quantity 
coincide. However, this fact does not necessarily hold for more general matter fields. 

For spherically symmetric homothetic spacetimes, we can assume that there is a coordinate system 
t and r such that all dimensionless variables are functions of a single dimensionless self-similar variable 
£ = r/t. The solution is invariant under scale transformation t = at, r = ar for any constant o. Thus 
the self-similar variables can be determined from dimensional considerations in the case of homothety. 
Therefore, we can conclude homothety as the general relativistic analogue of complete similarity. 

From the constraints (|3.9j) and (|3.1U|) . we can show that if we consider the barotropic equation of 
state, i.e., p = f(fJ,), then the equation of state must have the form p = K/j,, where K is a constant. 
This class of equations of state contains a dust fluid (K = 0), a radiation fluid (K = 1/3) and a stiff 
fluid {K = 1) as special cases. Other important matter fields that are compatible with homothety are 
a massless scalar field and a scalar field with an exponential potential. 
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3.2 Kinematic self-similarity 



Although homothetic solutions can contain several interesting matter fields, the matter fields compat- 
ible with homothety are rather limited. In more general situations, matter fields will have intrinsic 
dimensional constants. For example, when we consider a polytropic equation of state, such as p = K/i 7 , 
the constant K has dimension M 1- " 7 'L 3 ^ 7 " 1 ), where we should be reminded that we have chosen the 
light speed c to be unity. We can also consider a massive scalar field, where the mass of the scalar 
field has dimension M. In such cases, it is impossible to assume homothety because the system has 
a characteristic scale. By analogy, we can consider the general relativistic counterpart of incomplete 
similarity. From comparison with self-similarity for a polytropic gas in Newtonian gravity, kinematic 
self-similarity has been defined in the context of relativistic fluid mechanics as an example of incom- 
plete similarity |191 I2U1 I23j . It should be noted that the introduction of incomplete similarity to 
general relativity is not unique. For example, partial self-similarity has been defined and applied to 
inhomogeneous cosmological solutions jHZl EM 1^ • 

A spacetime is said to be kinematic self-similar if it admits a kinematic self-similar vector £ which 
satisfies the conditions 

Z^V = 2<5 V» ( 3 - u ) 
Cgu^ = au^, (3.12) 

where u M is the four-velocity of the fluid and = g^ u + u^u u is the projection tensor, and a and 
5 are constants |19l 1201 123 j . If 5 ^ 0, the similarity transformation is characterized by the scale- 
independent ratio a/S, which is referred to as the similarity index. If the ratio is unity, £ turns out 
to be a homothetic vector. In the context of kinematic self-similarity, homothety is referred to as 
self-similarity of the first kind. If a = and 5 ^ 0, it is referred to as self-similarity of the zeroth 
kind. If the ratio is not equal to zero or one, it is referred to as self-similarity of the second kind. If 
and 5 = 0, it is referred to as self-similarity of the infinite kind. If 5 = a = 0, £ turns out to be 
a Killing vector. 

From the Einstein equation (|3.5[) . we can derive 

C^G^ = SttGjC^. (3.13) 

This equation is called the integrability condition. Now we can rewrite the integrability conditions 
(|3.13j) in terms of kinematic quantities of the fluid. The covariant derivative of the fluid four velocity 
is decomposed into the following form: 

Ufj,-v = <V + -Oh^ + uj^ v - Ufj,u nu , (3.14) 

where 

V = VX)V;A, (3.15) 

6 = gTB^, (3.16) 







(3.17) 






(3.18) 


u 2 


1 


(3.19) 






(3.20) 
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where the semicolon denotes the covariant derivative. Using the above quantities, the integrability 
condition 1)3. 13j) is rewritten as follows (cf. p3 ): 



(5-a)(-8uj 2 



2n. 



8ttG 



1 3 

+ 2a//) + - (Cgp + lap) 



2(5 - a)(6 + 9 2 - Auj 2 ) = 8ttG 



hc^ + 25n)-hc ( :p + 25p) 



2uj X u,u^ + 2u 



4lu 2 u x = 0, 



(3.21) 

(3.22) 
(3.23) 



o\ p - u p (j\ u u u - uxfjppii 11 + 6(j\p + cr\ K LU K p + <T pK uJ K A + 2u A K u Kp + -h Xp uj 2 = 0. (3.24) 

For the first-kind case, in which a = 5 ^ 0, equations (|3.9jl and (|3.1Ujl are obtained from equations 
(|3.21|) and 1)3. 22|) . When a perfect fluid is irrotational, i.e., = 0, the Einstein equations and the 
integrability conditions (|3~2TI) - (|3~24ll give |23 E0] 



(a - 6)Kuv = 0, 



(3.25) 



where TZp U is the Ricci tensor on the hypersurface orthogonal to u^. This means that if a solution 
is kinematic self-similar but not homothetic and if the fluid is irrotational, then the hypersurface 
orthogonal to fluid flow is flat. 



4 Spherically symmetric self-similar solutions 

4.1 Spherically symmetric solutions 

Although self-similar solutions can play important roles even in nonspherically symmetric solutions, 
such as homogeneous cosmological models |4U1 173j , we focus in the rest of this article on spherically 
symmetric spacetimes. The line element in a spherically symmetric spacetime is given by 

where dQ, 2 = d9 2 + sirpQdtp 2 . We consider a perfect fluid as a matter field, for which the energy- 
momentum tensor is given by equation (|3.7|) . We adopt the comoving coordinates, where the four- 
velocity of the fluid u M has the components 

u„ = (-e*> 0,0,0). (4.2) 

Then, the Einstein equations and the equations of motion for the perfect fluid are reduced to the 
following simple form: 



(fj, + p)$ r = -p r , (4.3) 

(M + P)*t = -iH-2(ii + p)^, (4.4) 

m r = 4irfiR r R 2 , (4.5) 

m t = -4irpR t R 2 , (4.6) 

= -R tr + $ r Rt + V t Rr, (4.7) 

2Gm = R(l + e- 2 *R t 2 -e- 2 *R r 2 ), (4.8) 

where the subscripts t and r denote the partial derivatives with respect to t and r, respectively, and 
m(t, r) is called the Misner-Sharp mass. When a perfect fluid obeys an equation state p + fj, = 0, 
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which is equivalent to a cosmological constant, the first two equations are trivially satisfied. In this 
case, one can use the following equations: 



R 2 + 



Rt \ Rt T 
R R 



+ e 



2*-2<I' 



Rh Rt T / Rt 
R R \ R 



R R I R 



R r \ R r „ 



•>ii> / T o T T Rtt Rt^t Rt&t 



i? 



i? 



i? 



SirGpe 



2* 



8vrGpe 2 *, (4.10) 



< - 2 * f $ rr + $ r 2 - $ r ^ r + ^ + ^ - ^ ) - -swap. 



(4.c 



(4.11) 



which are (tt), (rr) and (00) components of the Einstein equations, respectively. Five of the above 
nine equations are independent. 



4.2 Spherically symmetric homothetic solutions 

There is a large variety of spherically symmetric homothetic solutions. The pioneering work in this 
area was done by Cahill and Taub The application contains primordial black holes ^Hl El El- 

cosmological voids H H3 [7UJ HQ , cosmic censorship jSHl EZl EH IZH ESI EOl HH HH1 EZI and critical 
behavior |21 | 1421 l25j . See [HE] and [HOI for recent reviews of cosmic censorship and critical behavior, 
respectively. The classification of all spherically symmetric homothetic solutions with a perfect fluid 
has been made (2H1 E2 El ESI ESI EH1 ■ The spacetime structure possible for homothetic solutions has 
been studied ^7j. The special case where the homothetic vector is orthogonal or parallel to the fluid 
flow has also been studied [23 122] • It has been revealed that a homothetic solution describes the 
dynamical properties of more general solutions in spherically symmetric gravitational collapse |37j . 
The stability of homothetic solutions has been studied mEimEHlEZlEHlHEE]- 

When the spacetime admits a homothetic vector, which is neither parallel nor orthogonal to the 
fluid flow, the homothetic vector £ can be written as 



£_ f d_ + r d_ 

dt dr ' 



and the self-similar variable £ is given by 



Homothety implies that the metric functions can be written 



(4.12) 



(4.13) 



(4.14) 



As we have seen, the equation of state must be of the form p = Kp for homothetic spacetimes. Then 
the governing equations for homothetic solutions are written as 



„2$ 



2$ 



M + M' 

M' 
M 
~~S 



4K 

a a i 1+K P 

2 



2K 
' 1+K 



pS 2 (S + S') , 
-KpS 2 S', 

1 + a' 1 (vC 2 ) fs' 2 - pw?S\S + S') 2 , 



(4.15) 
(4.16) 
(4.17) 
(4.18) 

(4.19) 
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where a a and a w are integration constants and the prime denotes the derivative with respect to ln£. 
The dimensionless functions r/(£) and M(£) are defined by 

SttG// = \ (4.20) 

2Gm = tM. (4.21) 

The above formulation is based on |111 15115], It is possible to choose another function in the same co- 
moving coordinates, as adopted in [1811131 HI] . In the comoving coordinates, the dynamical properties 
of the fluid elements are very clear. 

There are other useful formulations in analyzing homothetic solutions. One of the most natural 
coordinate systems for homothetic spacetimes is the so-called homothetic coordinates. In terms of 
this coordinate system, the dynamical systems theory has been applied to homothetic solutions with 
a perfect fluid for classification [71 1281 12*9*] . In the homothetic coordinates, the self-similar variable is 
chosen to be the spatial or time coordinate, depending on whether the homothetic vector is timelike 
or spacelike. If the homothetic vector is timelike, the line element is written as 

ds 2 = e 2t [-Dj(x)dt 2 + dx 2 + D%{x)dn 2 ] . (4.22) 

If the homothetic vector is spacelike, the line element is written as 

ds 2 = e 2x [-dt 2 + D 2 (t)dx 2 + D^(t)dQ 2 ] . (4.23) 

If the homothetic vector is timelike in one region and spacelike in another region of the same spacetime, 
the above two charts must be patched on the hypersurface on which the homothetic vector is null. 

Another coordinate system is that of area coordinates, in which the physical properties of the 
spacetime are clear. The area coordinate system has been adopted [57J EH1 EH) • In this coordinate 
system, the line element in homothetic spacetimes is written as 

ds 2 = -e 2 ^dt 2 + e 2 ^dr 2 + r 2 dn 2 , (4.24) 

(4.25) 



r 



t 

where u* and u r are also to be determined. 

4.3 Spherically symmetric kinematic self-similar solutions 

A kinematic self-similar vector may be parallel, orthogonal or tilted, i.e., neither parallel nor orthog- 
onal, to the fluid flow. Spherically symmetric kinematic self-similar perfect fluid solutions have been 
recently explored by several authors |21 131 ITU1 1351 150] . 

In a spherically symmetric spacetime, the kinematic self-similar vector field £ is written in general 

as 

t = hi{t,r)^ + h^r)^ (4.27) 

in the comoving coordinates, where h\(t,r) and hi{t,r) are functions of t and r. When hi = 0, £ is 
parallel to the fluid flow, while when h\ = 0, £ is orthogonal to the fluid flow. When both h\ and hi 
are nonzero, £ is tilted. 

When the kinematic self-similar vector £ is tilted to the fluid flow and not of the infinite kind, £ 
and the metric tensor g^ v are written in appropriate comoving coordinates as 

t = (at + 0)^ + ^, (4.28) 
ds 2 = -e 2 ^) dt 2 + e 2 ^) dr 2 + r 2 S(Z) 2 dQ 2 , (4.29) 
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where a is the index of self-similarity. For a = 1, i.e., homothety or self-similarity of the first kind, 
we can set (3 = and then £ is given by £ = r/t. For a = 0, i.e., self-similarity of the zeroth kind, we 
can set (3 = 1 and then £ is given by £ = r/e*. For a/0 and a / 1, i.e. self-similarity of the second 
kind, we can set (3 = and then £ is given by £ = r/(ai) 1 / a . If the kinematic self-similar vector £ is 
tilted to the fluid flow and of the infinite kind, £ and the metric tensor g^ u are written in appropriate 
comoving coordinates as 

« = 4t +r t (4 - 3o) 

ds 2 = _ e 2*(O dt 2 + ^_ dr 2 + (431) 

where the self-similar variable is given by £ = r/t. 

If the kinematic self-similar vector £ is parallel to the fluid flow and not of the infinite kind, we 
have in appropriate comoving coordinates as 

« = «|. 

A, 2 = -( 2 <°-V*<''><if 2 + fV + ( 2 S(r)W, (433) 

where a is the index of self-similarity and the self-similar variable is given by £ = r. If the kinematic 
self-similar vector £ is parallel to the fluid flow and of the infinite kind, we have in appropriate comoving 
coordinates as 

< = <|. <««> 

rfs 2 = -i^l'tftiJ+SH'af, (435) 

where the self-similar variable is given by £ = r. 

If the kinematic self-similar vector £ is orthogonal to the fluid flow and not of the infinite kind, we 
have in appropriate coordinates 

£ = r± (4.36) 
ds 2 = -r 2a dt 2 + e 2 ^dr 2 + r 2 S(t) 2 dn 2 , (4.37) 

where a is the index of self-similarity and the self-similar variable is given by £ = t. If the kinematic 
self-similar vector £ is orthogonal to the fluid flow and of the infinite kind, we have in appropriate 
coordinates 

£ = r± (4.38) 

P 2*(t) 

ds 2 = -r 2 dt 2 + —p^-dr 2 + S(tfd£l 2 , (4.39) 

where the self-similar variable is given by £ = i. 

Not as homothetic solutions in the tilted case, kinematic self-similar solutions have a characteristic 
structure. We now show an example of them in the case of self-similarity of the second kind, where a 
kinematic self-similar vector is tilted to the fluid flow. In this case, the Einstein equations imply that 
the quantities m, fi and p must be of the following form: 

1dm r 2 

" = M 1 (0 + ^M 2 (0, (4.40) 
SirGur 2 = W 1 (0 + ^W 2 (0, (4.41) 

2 

SirGpr 2 = P 1 (0 + ^P 2 (0, (4.42) 



r 

„2 
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where £ = r/(at) 1 ^ a . In other words, dimensionless quantities on the left hand side are decomposed 
into two parts, one remains constant and the other behaves as (r/t) 2 oc r 2<yl ~°^ as £ is fixed. Then, the 
original partial differential equations are satisfied when and only when the Einstein equations and the 
equations of motion for the matter field are satisfied for each of the O(l) and 0[(r/t) 2 ] terms. The 
equations ()4.3|) - (|4.1Uj) for a perfect fluid then reduce to the following: 



M 1 +M[ = W^iS + S'), (4.43) 

2>M 2 +M' 2 = W 2 S 2 (S + S'), (4.44) 

M[ = -Pi5 2 5', (4.45) 

2aM 2 +M 2 = -P 2 S 2 S\ (4.46) 

Mi = S[l - e~ 2 * {S + S') 2 ], (4.47) 

a 2 M 2 = SS'V 2 *, (4.48) 

(Pi + Wi)*' = 2P 1 -P[, (4.49) 

{P 2 + W 2 )& = -P' 2 , (4.50) 

W{S = -(Pi + W 1 ){^'S + 2S'), (4.51) 

(2aW 2 + W 2 [)S = -(P 2 + W 2 )(*'S + 2S'), (4.52) 

S" + S' = S'& + (S + 5')*', (4.53) 

S'(S' + 2^'S) = a 2 W 2 S 2 e 2 ^, (4.54) 

2S(S" + 2S')-2^'S{S + S') = -S' 2 - S 2 + e 2 *(l-W 1 S 2 ), (4.55) 

2S(S" + aS' -^'S') + S' 2 = -a 2 P 2 S 2 e 2 ^, (4.56) 

(S + S'){S + S' + 2&S) = (l + PiS 2 )e 2 *, (4.57) 



where we have omitted the bars of $ and ^ in (|4.29|) for simplicity and the prime denotes the derivative 
with respect to ln£. A similar structure of basic equations can be found for kinematic self-similar 
solutions of the second, zeroth and infinite kinds both in the tilted and orthogonal cases and of the 
second and zeroth kind in the parallel case. The exceptions are the first kind in the tilted, parallel 
and orthogonal cases and the infinite kind in the parallel case. See |49( 150) for the basic equations for 
spherically symmetric self-similar solutions for all cases. 

It is interesting to consider the spherically symmetric self-similar solutions of the infinite kind 
with a kinematic self-similar vector parallel to the fluid flow. The metric form demanded by this 
self-similarity, which is given by equation ()4.35j) . is nothing but the general form of the line element 
in spherically symmetric static spacetimes when the chosen radial coordinate is the radial physical 
length. Therefore, all static solutions have a kinematic self-similar vector of the infinite kind that is 
parallel to the fluid flow. Inversely, all spherically symmetric solutions with a kinematic self-similar 
vector of the infinite kind parallel to the fluid flow are static. The equation of state is not restricted 
at all. 

4.4 Equation of state 

It is obvious that equations (|4.41j) and (|4.42|) strongly restrict the form of the possible equations of 
state. The detailed analysis shows the following restriction on the equation of state [35]. Suppose we 
have the barotropic equation of state, i.e., p = /(/it). The self-similarity of the second kind with the 
index a can admit only the following equation of state: 

k lX + k 2 x a = f(C lX + C 2 x a ), (4.58) 

where k±, k 2 , C\ and C 2 are arbitrary constants. For the self-similarity of the zeroth and infinite kinds, 
we cannot determine the equation of state alone from the decomposed form. It should be noted that 
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the above discussion is based not on the whole equations but only on the decomposed form of p and 
p such as equations (|4.41|) and (|4.42|) . 

For later convenience, we introduce and focus on the following equations of state: 

• Equation of state (1) (EOS1) 

p = Kff, (4.59) 
where K and 7 are constants. Here we assume that K 7^ and 7 7^ 0, 1, 

• Equation of state (2) (EOS2) 

{p = Kn 1 , 
. P (4.60) 
p = m b n-\ -, ^ > 

7-I 

where the constant and n(t, r) correspond to the mean baryon mass and the baryon number 
density, respectively. Here we assume that K 7^ and 7 ^ 0, 1. In the literature, this equation 
of state is sometimes called a relativistic polytrope, 

• Equation of state (3) (EOS3) 

p = Kp. (4.61) 

where we assume that — 1 < K < 1. 

EOS1 and EOS2 are two kinds of polytropic equations of state. These equations of state are incom- 
patible with homothety. For < 7 < 1, both EOS1 and EOS2 are approximated by a dust fluid in the 
high-density regime. For 1 < 7, EOS2 is approximated by EOS3 with K = 7 — 1 in the high-density 
regime. For 2 < 7 for EOS2 and 1 < 7 for EOS1, the dominant energy condition can be violated in 
the high-density regime, which would be unphysical. 



5 Exact spherically symmetric self-similar solutions 
5.1 Vacuum 

In a vacuum, the only spherically symmetric solutions are the Minkowski solution and the Schwarzschild 
solution from Birkhoff 's theorem. Both solutions have kinematic self-similar vectors. Although there 
are no fluids, we can introduce a unit timelike vector u^. 

The Minkowski solution has seven kinematic self-similar vectors including a homothetic vector in 
the tilted case. The Minkowski solution is represented by 

ds 2 = -dt 2 + dr 2 + r 2 dQ 2 , (5.1) 
2Gm = 0, (5.2) 
8irGp = 8irGfi = 0. (5.3) 

The metric can be represented in the Milne form 

ds 2 = -dr 2 + r 2 dp 2 + t 2 sinh 2 pdfi 2 , (5.4) 

where t = r cosh p and r = r sinh p, or in another form 

ds 2 = —m 2 dv 2 + dm 2 + w 2 cosh 2 vdtt 2 , (5-5) 

where t = wsinhi/ and r = tucoshzA This spacetime has the following kinematic self-similar vectors: 
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• First kind, tilted 



• First kind, parallel 



• First kind, orthogonal 



Second kind with any a, tilted 



Zeroth kind, tilted 



Zeroth kind, orthogonal 



• Infinite kind, parallel 



d d 



d 



ZD 



_d_ 

dm' 



at— r— 
dt <9r ' 



8_ d_ 

dt dr ' 



dr ' 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



The Schwarzschild solution has two kinematic self-similar vectors but does not have a homothetic 
vector. The Schwarzschild solution is represented by 

,2 

(5.13) 



ds 2 = -fl-^W + - ^ 



2GW 



+ rW, 



2Gm = 2Gm , 
8irGp = 8nGp = 0, 



(5.14) 
(5.15) 



where uiq is a constant. This spacetime can be represented in the following choice of coordinates: 

dp 2 



ds 2 = -dr 2 + (2Gm )^i- 

\[Up- t )\ 

This spacetime has the following kinematic self-similar vectors: 







4/3 \ 


+ 




dn 2 



(5.16) 
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Second kind with a = 3/2, tilted 



T It +P W (5 - 17) 



Infinite kind, parallel 



4- (5.18) 

5.2 Cosmological constant 

Since the cosmological constant introduces a length scale l/-y/|A|, solutions cannot be homothetic. 
However, the de Sitter solution, the Schwarzschild-de Sitter solution and the Nariai solution admit 
kinematic self-similar vectors. 

The de Sitter solution is represented by 

ds 2 = -dt 2 + e 2 v^/3*(dr 2 + rW), (5.19) 

2Gm = ^rVv 7 ^*, ( 5 . 20 ) 

8ttGp = SirGp = -8ttGA, (5.21) 

where A is a cosmological constant. This solution is represented in the static coordinates as 

ds 2 = - (l - l^dr 2 + (l - ^\/) \p 2 + p 2 dn 2 , (5.22) 

2Gm = -Ap 3 , (5.23) 

8vrGp = SttGp = SttGA. (5.24) 

When A is negative, the solution is called the anti de Sitter solution. The de Sitter solution has the 
following kinematic self-similar vectors: 

• Zeroth kind, tilted 

¥t + x 4- (5 ' 25) 

where A is a non-zero constant. 

• Zeroth kind, parallel 

!• < 5 - 26 > 

• Zeroth kind, orthogonal 

r± (5.27) 
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Infinite kind, parallel 



dr' 



(5.28) 



The Schwarzschild-de Sitter solution is an exact solution with a cosmological constant, which is 
represented by 



ds 2 

2Gm 
8irGp 



1 



-i 



2Gm ° 1 \p 2 )dr 2 I (l " A7 "'" l ' :2] •' • 2 1 :2 '"' 2 



2Gm + -Ap 3 , 
-8ttG/x = -8vrGA. 



2(7 '"" -Ap 2 ) d^ + pW, (5.29) 
P 3 / 

(5.30) 
(5.31) 



where tuq is a constant and A is a cosmological constant. When A is negative, the solution is called the 
Schwarzschild-anti de Sitter solution. The Schwarzschild-de Sitter solution has the following kinematic 
self-similar vector: 

• Infinite kind, parallel 



dr 



(5.32) 



The Nariai solution is an exact solution with a cosmological constant, which is represented by 

2 „ 1 



ds 2 = - a(t) sin (ln(VXr)J + b(t) cos (ln(VXr)J dt 2 + j^idr 2 + r 2 dQ?) 

2Gm = 1/VA, 
8vrGp = SirGp = -A, 

where a and 6 are arbitrary functions of t. With the choice 

t i/( Cl VX)-i) 



ciVX 



■(^cos(lnt) + J Bsin(lnt)), 



(-^sin(lnt) + Scos(lnt)), 
ciVA 

where A and .B are constants, and the coordinate transformation 



(5.33) 

(5.34) 
(5.35) 



(5.36) 
(5.37) 



this metric is written as 

ds 2 = _ J_i2/( Clv ^)-2 

cfA 



A 



A sin I In — | + B cos (In — 
t V t 



dt 2 + -^(dr 2 + r 2 dn z ). 
Ar 2 



This spacetime is also written in the static coordinates as 



ds' 



Asm (ln(\/Ar) ) + B cos (ln(VXr) 



dr2+ A^ (dr2 + rW) ' 



(5.38) 



(5.39) 



(5.40) 



where A and B are constants. This spacetime has the following kinematic self-similar vectors: 
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• Infinite kind, tilted 



'i*'*?- (541) 



Infinite kind, parallel 



dr 



(5.42) 



5.3 Dust fluid 

Without assumption of self-similarity, the general solution for a spherically symmetric dust fluid is 
exactly obtained, which is called the Lemaitre-Tolman-Bondi (LTB) solution [461 1661 15]. 

Therefore, spherically symmetric self-similar solutions with a dust fluid are a subclass of the LTB 
solutions. Homothetic solutions in the tilted case was completely classified j!2j . The homothetic LTB 
solutions have been discussed in the context of cosmic censorship jSH]. The homothetic LTB solution 
is represented by 



ds 2 
2Gm 



.<ft2 + ^ + ^f dr 2 + r 2 S 2 dQ 2 , 
1 + 2E 



2V1 + 2Er, 



2£r 



2 S 2 (£S± ^2E + 2T/S)' 



(5.43) 
(5.44) 
(5.45) 



where £ = r/t, where a subscript £ means the derivative with respect to £, and where E and T are 
constants with a relation 



E 



1 



(r 2 



i) 



(5.46) 



S is given by 



(2E)\ 



VES 2 + TS 2T_ 

3 6 

2r 

(-2E)\ 



■ sinh 



ES 

r 



for E > 0, 



■ sin 



-ES VES 2 + TS 



V2E 



where D is a constant. This solution has a homothetic vector 

d d 
t— + r— . 
at or 



for E 
for — 



0, 



(5.47) 



\ < E < 0, 



(5.48) 



E can be interpreted as the sum of the kinetic energy and the potential energy per unit mass. When 
E = 0, each shell is marginally bound. This solution is a two-parameter family of solutions of E and 
D and reduces to the flat Friedmann-Robertson- Walker (FRW) solution when E = D = 0. It is noted 
that there are nonmarginally bound LTB solutions with a homothetic vector. 

In addition, there are kinematic self-similar solutions of the second, zeroth and infinite kinds. From 
equation Q3.25JI . any three-surface orthogonal to the fluid flow in a kinematic self-similar solution to 
the Einstein equations of the second, zeroth or infinite kind that contains only irrotational dust as a 
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matter field is flat. Because the flatness of the three-surface implies that these solutions are marginally 
bound, a spherically symmetric kinematic self-similar solution to the Einstein equations of the second, 
zeroth or infinite kind that contains only a dust fluid is described by the marginally bound LTB 
solutions. These solutions have been investigated by several authors [T9l f2| fTTTl lB0| . 
The second-kind kinematic self-similar LTB solution is represented by 

ds 2 - dt 2 I 9[^/(3-^) + (2a/3-l)^ 2 r2 , 3a /(3~ 2a ) , |4 /3 d o 2 (54Q) 

87TG V = 9 [ Kr 3a/(3-2 Q ) + (2a ( /3 _ 1) ) t](Kr 3 a /(3-2 a) _ f) > ( 5 ' 50 ) 

4 o 

2Gm = -r 3 , (5.51) 
9 

where k is an arbitrary dimensional constant and a ^ 3/2. For a = 3/2, the solution turns out to be 
the flat FRW solution. This spacetime has the following kinematic self-similar vector: 

8 3 — 2a d 

at m + — r lTr- (5 ' 52) 

The zeroth-kind kinematic self-similar LTB solution is represented by 

ds 2 = -dt 2 + {t ~ 2A/ N 3 ~ V + r 2 \t - Alnr| 4 /W, (5.53) 
\t — Xmr] 2 ' 6 

^ = 3(t-Alnr)(t-2A/3-Alnr) ' (5 ' 54) 

2Gm = ^r 3 . (5.55) 

This spacetime has the following kinematic self-similar vector: 

where A is an arbitrary dimensional constant. 

The infinite kind kinematic self-similar LTB solution is represented by 

ds 2 = -dt 2 + t 2 \ar~ 3 / 2 -t\~ 2/3 dr 2 + r 2 \ar~ 3 / 2 -t^dft 2 , (5.57) 

2Gm = ^r 3 . (5.59) 

This spacetime has the following kinematic self-similar vector: 

8 2 8 

where a is an arbitrary dimensional constant. 

The flat FRW solution with a dust fluid will be discussed together with those with a perfect fluid. 

5.4 Perfect fluid 

5.4.1 Homothetic solutions 

As examples of homothetic exact solutions with a perfect fluid obeying the equation of state p = Kfj,, 
we discuss the power-law flat FRW solution, the homothetic static perfect fluid solution, and the 
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Kantowski-Sachs solution. The power-law flat FRW solution and the homothetic static solution have 
kinematic self-similar vectors as well as a homothetic vector. The Kantowski-Sachs solution has no 
kinematic self-similar vector except a homothetic vector. 
The flat FRW solution has the following form: 

ds 2 = _ rft 2 + ^4/(3(1+^)] (dr 2 +r 2 dn2)j (5 61) 

2Gm = (5.62) 

AK 

8nGp = *«GKp= 3{l + K)2t2 . (5.63) 

This solution has the following kinematic self-similar vectors: 
• First kind, tilted (K / -1/3) 

d 1 + 3K 8 
t^ + o„ ■ *rS —, (5.64) 



dt 3(1+ if) dr 

First kind, parallel (K = -1/3) 



dt 

Second kind for a / 3(1 + K)/2, tilted 



4, (5-65) 



'dt V 3(1+ if)/ dr' 
• Second kind with a = 3(1 + K)/2, parallel 



«4 + f 1 - 7777^7^ r#, (5.66) 



• Zeroth kind, orthogonal 



aij^, (5.67) 



r|-, (5.68) 



Infinite kind, tilted 



4t - WTKj^r- (5 ' 69) 
The homothetic static perfect fluid solution is represented as the following: 

ds 2 = -r 4K '< 1+K >dt 2 + K2 ,^kV + rW ' (5 ' 70) 

2Gm = 7OTTT r ' 2 (5T1) 

8. G p = 8,G^ = (g2+ g +1)r2 . (5.72) 

Since the center r = is singular and timelike, it must be a naked singularity. This solution has the 
following kinematic self-similar vectors: 
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• First kind, tilted (K / 1) 



1 + K dt dr 

First kind, orthogonal (K = 1) 



1-K d d . „ , 

^ + ^- (5-73) 



r±, (5.74) 



Second kind for a / 2if/(l + K), tilted 



2K \ ^ d 



Second kind with a = 2if/ (l + K), orthogonal 



r|-, (5.76) 



• Zeroth kind, tilted 



2 ^ ^ d 

*m ( 5 - 77 ) 



l + K St <9r 
The Kantowski-Sachs solution is represented by 

^ 2 = -eft 2 + f^/^dr 2 + (1 ^ff. 1} ^ 2 > (5-78) 
4(l + K)K 2 t 

2Gm = V 'T / — , (5.79) 



(1 + 3K) 3 / 2 (K- l) 3 / 2 ' 
4if 2 
(1 + K) 2 t 2 ' 



4if 2 

8ttGp = 8irGK[i = — (5.80) 



with a homothetic vector 



^ 1 + 3K d , rtvl , 

'S + TTx'S- (5 ' 8l) 

— 1 < if < —1/3 must be satisfied for this solution to be physical and in that case the above homothetic 
vector is tilted. 



5.4.2 Nonhomothetic kinematic self-similar solutions 

As examples of nonhomothetic kinematic self-similar solutions, we discuss the general FRW solutions 
and the Gutman-Bespal'ko solution. 
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The general FRW solution has kinematic self-similar vectors and does not have homothetic vectors 
if it is not power-law or if it is not flat. The general FRW solutions are given by 



ds 2 = -dt 2 + a(t) 2 {dr 2 + S{r) 2 d£l 2 ), (5.82) 
2Gm = aS(l -S' 2 + a 2 S 2 ), (5.83) 

8ttG P = -2^-(^V-4, (5.84) 




8itGh = 3[-) + — , (5.85) 

for k = 1 

S(r) = { r, for k = (5.86) 

for k = —1 

where a dot and a prime denote the derivatives with respect to t and r, respectively. 

The non-power-law flat FRW solution (k = 0) has the following kinematic self-similar vector 
independent of the form of the equation of state: 

• Zeroth kind, orthogonal 

r|-. (5.87) 

The closed FRW solution (k = 1) has the following kinematic self-similar vector for EOS2 with 
7 = 2/3: 

• Second kind with a = 3/2, parallel 

4- (5 - 88) 

The curved FRW solutions {k = ±1) with the equation of state p = — fi/3 have the following 
additional kinematic self-similar vector: 

• First kind, parallel 

•s +r l- (5 - 89) 

Only for a stiff fluid p = /i, the Gutman-Bespal'ko solution exists, which is represented by |31 [ I64 | 

ism 

ds 2 = -^r 2 dt 2 + dr 2 + ^r 2 (l + a 1 e t + a 2 e~ t )dn 2 , (5.90) 

IGm = -^r=(l - 4oio 2 )(l + aie l + a 2 e~ t )~ 3/2 r, (5.91) 
2y 2 

8irGp = 8vrG^ = (l-4aia 2 )(l + aie* + a2e"*)" 2 r- 2 , (5.92) 

where a\ and a 2 are arbitrary constants. In this spacetime the physical center r = is singular. This 
spacetime has the following kinematic self-similar vector: 
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• First kind, orthogonal 



r-. (5.93) 

This solution includes the homothetic static perfect fluid solution for a stiff fluid as a special case 
o-i = a 2 = 0. 

The results in this section are summarized in Tabled As we will see in section H3 it can be shown 
that this table provides the complete list of kinematic self-similar solutions compatible with EOS1, 
EOS2 and EOS3. 



6 Nonexistence of kinematic self-similar solutions with a polytropic 
equation of state 

Among kinematic self-similarities, homothety in the tilted case includes a large variety of solutions 
and has been intensively investigated. As we have seen, the equation of state is restricted to be of the 
form p = Kfi for homothetic solutions. 

In this section, based on |481 1491 15U] . we will briefly see that there are no kinematic self-similar so- 
lutions with nontrivial polytropic equations of state although we first expected that the generalization 
from homothety to kinematic self-similarity would enable us to analyze a wider class of physical solu- 
tions. As we have already mentioned, it is possible to construct a characteristic length scale from given 
dimensional constants in the general relativistic system of a perfect fluid with a polytropic equation 
of state. According to a usual procedure, we have introduced incomplete similarity into this system. 
Kinematic self-similarity is a natural generalization of incomplete similarity into general relativity. 
Therefore, it is highly nontrivial that there are no kinematic self-similar solutions with a kinematic 
self-similar vector tilted to the flow of a perfect fluid with a polytropic equation of state. 



6.1 Tilted cases 
6.1.1 Second kind 

We consider the second-kind kinematic self-similar solutions in the tilted case, in which the energy 
density and the pressure of a perfect fluid are of the form (|4.41|) and (|4,42|) . respectively. The equation 
of state gives the relation among the functions Pi, P 2 , W\ and W 2 . EOS1 admits the following two 
cases: 

a = 7, P 1 = W 2 = 0, P 2 = {87TG p-i^ ~ 2lW l ( A ) 

P2 = W l= 0, Pl = ^.^ eWl (B) (6.2) 
while EOS2 admits the following two cases: 

a = 7 , P 1 = 0, P 2 = K 2 r 27 WT = (7 ~ 1)W 2 , (C) (6.3) 

« = -, ^2 = 0, P x = * = (7 - l)Wi. (D) (6.4) 

7 m b '(87rG)T 1 7^ 

We can show that none of these cases satisfies the Einstein equations although they are compatible 
with the decomposition given by equations (|4.41|) and (|4.42[) . Subtracting equation (J4.57I) from (|4.55f) 
and eliminating S" by use of equation (|4.53|) . we obtain 

2& = {P 1 + W l )e™ . (6.5) 
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Then equations (|4~i9")) and (|4~5U)) result in 



e 



2 *(Pi + I^i) 2 = 4P 1 -2P{, (6.6) 
e^(P l + W l ){P 2 + W 2 ) = -1P' 2 . (6.7) 

It is obvious that P\ = implies W\ = from equation (|6,6j) while P 2 = implies W2 = or 
Pi + Wi = from equation ()6.7|) , Therefore, it is concluded that all cases (A)-(D) result in vacuum 
spacetimes. 

6.1.2 Zeroth kind 

Next, we consider the zeroth-kind kinematic self-similar solutions in the tilted case. In this case, the 
Einstein equations imply that the quantities fi,p and m are of the forms 

9 Cm 

— = M 1 (£)+r 2 M 2 (0, (6.8) 
r 

8irGfir 2 = Wi(0 + r 2 W 2 (£), (6.9) 

SnGpr 2 = Pi(C) +r 2 P 2 (0, (6.10) 

where £ = r/e . A set of ordinary differential equations is obtained when it is stipulated that the 
Einstein equations and the equations of motion for the matter field be satisfied for the 0(1) and 0(r 2 ) 
terms separately. See |49| 150] for the complete set of the ordinary differential equations. Both from 
EOS1 and EOS2 with equations (jol)l) and lf6~TU|) . Pi = W\ = is concluded such that 

P 1 = Wi = 0, P 2 = K^itGf^Wi 1 , (A) (6.11) 

for EOS1, while 

P 1 = Wl = 0, P 2 = K 7 (V 2 --^-V, (B) (6.12) 
(87rG) 7 i m 6 ' \ 7 — 1 J 

for EOS2. From equations (|4.5]1 . (|4.6j) and (|4.8|) . we obtain 

e 2 * = (S + 5') 2 , (6.13) 
3,S,S'V 2 * = P 2 S 2 S' + W 2 S 2 {S + S'), (6.14) 

where we have omitted the bars of and in (]4.29|) for simplicity and the prime denotes the derivative 
with respect to ln£. Equation (|4,10|) gives 

2S(S" - &S') + S' 2 = -P 2 5 2 e 2 *, (6.15) 
{S + S')(S + S' + 2&S) = e 2 *. (6.16) 

exp($) = Co is concluded from equations (|6,16|) and (|6,13|) . where cq is a positive constant. Then 
P 2 = po is obtained from equation (]4.3|) . where po is a constant, which implies that W 2 = wq, where 
wq is a constant. Then, equation (|6.14|) gives the evolution equation for S: 

^(%) -(Po + ^o)%-^o = 0. (6.17) 



S ' S 



The solution to this equation is S = so£ 9 , where sq and q are constants, q / -1 must be satisfied 
because of equation (|6.13|) . Equation (|4.4f) with the fact P 2 + W 2 ^ gives 

^'S + 2S' = 0. (6.18) 

Then, the equality q = can be obtained from equation (|6.13|) and (|6,18|) . which implies that S = so- 
Finally, equations (|6.15|) and (|6.17|) give po = and wo = 0, respectively. Therefore, it is concluded 
that cases (A) and (B) result in vacuum spacetimes. 
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6.1.3 Infinite kind 



Finally, we consider the infinite-kind kinematic self-similar solutions in the tilted case. In this case, 
the Einstein equations imply that the quantities fj,,p and m are of the forms 

2Gm = Mi(0A 2 + M 2 (^), (6.19) 
8irGfi = Wi(£)/i 2 + W2(£), (6.20) 
SirGp = Pi(0/t 2 + P2(0, (6-21) 

where £ = r/t. A set of ordinary differential equations is obtained when it is demanded that the 
Einstein equations and the equations of motion for the matter field be satisfied for the O(l) and 
0(t~ 2 ) terms separately. See |49| 150) for the complete set of ordinary differential equations. Both 
from EOS1 and EOS2 with equations (l6~2U]l and (KT2T|) . P x = W x = is concluded such that 

P x = Wx = 0, P 2 = K^G^W^, (A) (6.22) 

for EOS1, while 

Pl = Wl = 0, P 2 = K 7 (V 2 --^-V, (B) (6.23) 

for EOS2. From equations (|4T7J) . (|4*U)) and (|4~TU)) . we obtain 

5" = + (6-24) 

(l-W 2 S 2 )e 2 * = 2SS" + S' 2 - 2^'S'S, (6.25) 
(1 + P 2 5 2 )e 2 * = S'(2$'S + S'), (6.26) 

where we have omitted the bars of $ and \D r in (j4.31j) for simplicity the prime denotes the derivative 
with respect to ln£. From equations (|6.24j) . (|6.25j) and ()6.26|) . 

P 2 + W 2 = 0, (6.27) 

is obtained, which implies that p = —fi and gives a contradiction. Therefore, it is concluded that cases 
(A) and (B) result in vacuum spacetimes. 

As a result, it is shown that there is no kinematic self-similar solutions with a nontrivial polytropic 
equation of state in the tilted case. 



6.2 Nontilted cases 

Even when the parallel and orthogonal cases are considered, except for the infinite-kind kinematic 
self-similar solutions in the parallel case which include all static solutions, the only possible solutions 
are the flat FRW solution as a zeroth-kind kinematic self-similar solution in the orthogonal case both 
for EOS1 and EOS2 and the closed FRW solution as a second-kind kinematic self-similar solution with 
an index a = 3/2 in the orthogonal case for EOS2 with 7 = 1/a = 2/3. 



7 Summary 

Self-similarity has been applied to many aspects of physics and other scientific fields. The introduction 
of self-similarity into Newtonian gravity is straightforward because it postulates absolute space and 
time. Since Newtonian gravity has only one dimensional constant, i.e. the gravitational constant, 
we can incorporate a polytropic gas as well as an isothermal gas into the framework of complete 
similarity. The introduction of self-similarity into general relativity is, however, not so straightforward 
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because there is no preferred coordinate system in the theory. The covariant definition of complete 
similarity in general relativity is homothety. Moreover, since two dimensional physical constants, the 
gravitational constant and the speed of light, are included in general relativity, it is impossible to 
incorporate many physically interesting matter fields, such as a polytropic equation of state, into 
the framework of homothety. This naturally leads to the introduction of incomplete similarity in 
general relativity. One of the most natural definitions of incomplete similarity in the fluid system in 
general relativity is kinematic self-similarity. Many known exact solutions turn out to be kinematic 
self-similar. At first glance it seems possible to construct kinematic self-similar solutions with a 
polytropic equation of state. However, more comprehensive study of the Einstein equations reveals 
that there are no such solutions. Although the present discussion implies somewhat limited application 
of kinematic self-similarity, there still remains a large possibility that kinematic self-similar solutions 
describe interesting gravitational phenomena of physically important matter fields, such as a double 
fluid system [2*3] , 
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Table 1: The complete list of kinematic self-similar solutions with a perfect fluid in the spherically 
symmetric spacetime for a perfect fluid with EOS1, EOS2 and EOS3. It is assumed that the energy 
density of the perfect fluid is not negative. See text and references therein. 



Matter field 


Kind 


Solution 


Vacuum 


first, parallel 


Minkowski 




first, orthogonal 


Minkowski 




second, tilted, any a 


Minkowski 




second, tilted, a = 3/2 


Schwarzschild 




zeroth, tilted 


Minkowski 




zeroth, orthogonal 


Minkowski 




infinite, parallel 


Minkowski 
Schwarzschild 


Cosmological constant 


zeroth, tilted 


de Sitter 




zeroth, parallel 


de Sitter 




zeroth, orthogonal 


de Sitter 




infinite, tilted 


Nariai 




infinite, parallel 


de Sitter 
Nariai 

Schwarzschild-de Sitter 


Dust 


first, tilted 


see Ref. Q2] 




second, tilted 


KSSLTB 




second, parallel, a = 3/2 


Flat FRW 




zeroth, tilted 


KSS LTB 




zeroth, orthogonal 


Flat FRW 




infinite, tilted 


KSS LTB 


Perfect fluid: EOS1 


zeroth, orthogonal 


Flat FRW 




infinite, parallel 


All static solutions 


Perfect fluid: EOS2 


second, parallel, a = 3/2 


Closed FRW with 7 = 2/3 




zeroth, orthogonal 


Flat FRW 




infinite, parallel 


All static solutions 


Perfect fluid: EOS3 


first, tilted 


see Refs. HE3EE31 




first, parallel 


FRW (K = -1/3) 




first, orthogonal 


Gutman-Bespal'ko {K = 1) 




second, tilted, a / 3(1 + K)/2 


Flat FRW 




second, tilted, a j= 2K/(1 + K) 


Homothetic static 




second, parallel, a = 3(1 + K)/2 


Flat FRW 




second, orthogonal, a = 2K/(1 + K) 


Homothetic static 




zeroth, orthogonal, 


Flat FRW 




infinite, tilted, 


Flat FRW 




infinite, parallel 


All static solutions 
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